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Approximate asymptotic conditions on the motion of compact, electrically charged particles are de- 
rived within the framework of general relativity using the Einstein, Infeld, Hoffmann (EIH) surface 
integral method. While superficially similar to the Abraham-Lorentz and Lorentz-Dirac equations, 
these conditions differ from them in several fundamental ways. They are not equations of motion in 
the usual sense but rather a set of conditions which the these motions must obey asymptotically in 
the future of an initial starting time. And furthermore, they do not admit the run-away solutions 
of these other equations. As in the original EIH work, they are integrability conditions gotton from 
integrating the empty-space (i.e., sourceless) Einstein-Maxwell equations of general relativity over 
closed two-surfaces surrounding the sources of the fields appearing in these equations.. No additional 
ad hoc assumptions, such as the form of a force law or the introduction of inertial reaction terms, 
are required for this purpose. Nor is there a need for any infinite mass renormalizations such as 
are required in other derivations since all integrals are over surfaces and thus finite. In addition to 
being asymptotic, the conditions of motion derived here are also approximate and apply, as do the 
original EIH equations, only to slowly moving systems. A 'slowness' parameter e is identified as the 
ratio of the light travel time across the system divided by a characteristic time, e.g., a period. Use is 
made of both the method of matched asymptotic expansions and the method of multiple time scales 
to obtain an asymptotic expansion in e and the expansion is carried to sufficiently high order to 
obtain the lowest order radiation reaction terms. The resulting conditions of motion are shown to 
not allow run-away motions. 

PACS numbers: 03. SOD, 04.40, 04.20M 



I. INTRODUCTION 

One of the most profound and far reaching consequences of Einstein's general theory of relativity is the restrictions 
it places on the motion of the sources of the fields described by the theory. This is a unique feature of general 
relativity and is found in no other field theory. In electrodynamics, for example, one must postulate not only the form 
of the Lorentz force law but also the form of the inertial forces that appear in the equations of motion. Furthermore, 
these restrictions, which I have called conditions of motion, follow from the empty-space field equations of the theory 
without the need to introduce source terms into these equations. As a consequence, many of the well-known difficulties 
and inconsistencies associated with point source models are avoided in general relativity. And finally, these conditions 
of motion are sufficiently restrictive that they limit the possible trajectories of an N-body system of compact sources 
to a 6N-parameter family of curves. In the final analysis, the origin of these conditions can be traced back to the 
general covariance of the field equations of general relativity, jl] 

Einstein and his co-workers Infeld and Hoffmann (EIH) first derived approximate conditions of motion for systems 
of slowly moving {v/c <^ 1) compact gravitating sources up to what is now called the post-Newtonian approximation 
in 1939. Q Two subsequent papers on the subject were pubfished by Einstein and Infeld in 1939 and 1949. [| The 
procedure used by EIH involved the evaluation of integrals of the field equations over surfaces surrounding the sources 
of the fields appearing in these equations. By using only surface integrals they were thus able to avoid the need to 
perform volume integrals over regions containing the sources. Even though EIH spoke of "point source" solutions 
it is clear that they were really refering to the form of the solutions in the region exterior to these sources so that 
their derivation was applicable to compact sources and black holes. While the EIH surface integrals were exact their 
evaluation required what has come to be called the slow motion approximation. In this approximation one assumes 
that the fields and source variables depend on time scales that are long compared to the light travel time across the 
system of sources. In a 1986 paper |^ (referred to hereafter as AI) I extended this idea by introducing multiple time 
scales into the theory and was thereby able to avoid the introduction of the fictitious dipole fields used by EIH. I also 
showed how one could make use of the method of matched asymptotic expansions to include the effects of radiation 
reaction into the conditions of motion. These conditions followed as integrability conditions directly from the empty- 
space Einstein field equations and were applicable both to sources with and without strong internal gravity and black 
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holes. Unfortunately, the EIH papers are referenced only in passing even in standard texts on general relativity |5[| 
and hardly at all in the literature. And yet these papers contain what I would claim is one of Einstein's greatest 
contributions to physics and the only reliable method to date for deriving conditions of motion. 

In this paper I have applied the EIH formalism together with the methods developed in AI to the case of electrically 
charged sources up to and including the order in which radiation reaction terms first appears. Past attempts to derive 
an expression for the electrodynamic radiation reaction force within the framework of Newtonian or special rclativistic 
physics have been beset with a number of well-known difficulties. The best known of the equations, the Abraham- 
Lorentz (AL) equation and its special rclativistic generalization, the Lorentz-Dirac (LD) equation (collectively the 
ADL equations) have been the source of continued discussion in the scientific literature. [|j - . All derivations of 
these equations have required one or more ad hoc assumptions concerning the form of a force law, an inertial reaction 
force or an indeterminate integral. Many derivations also require some form of infinite mass renormalization and 
they all make an arbitrary choice of the retarded over the advanced potential in constructing the field of a moving 
charge. Thus Dirac was forced to assume the form of the indeterminant (in fact, infinite) integrals that appeared in 
his derivation. 

But perhaps the most serious problem associated with the ADL equations is the existence of the famous (or 
infamous, depending on one's point of view) run-away solutions possessed by these equations. Their existence alone 
should have been sufficient evidence that there was something basically wrong with the whole enterprise. Dirac 
sought to eliminate these solutions by imposing the condition that in the asymptotic future all accelerations should 
be finite. This condition however leads to the phenomena of "pre-acceleration" with its attendant acausality. As 
a consequence, several authors have been led to propose alternative equations of motion for radiating charges. |l^] 
Such ad hoc modifications, motivated by little more than the desire to avoid the run-away sloution problem can 
hardly be considered to be a satisfactory solution to the problem. But more, it is far from obvious nor has it been 
shown that any of these equations are compatible with the restrictions imposed by general relativity on the motion 
of sources of the electromagnetic field. And indeed, if one accepts general relativity as a fundamental theory, one has 
no alternative than to accept, as one of its consequences, the restrictions it imposes on the motion of sources of both 
the electromagnetic and gravitational fields. 

The first derivation of conditions on the motion of electrically charged sources that did not involve an infinite mass 
renormalization and/or ad hoc assumptions of the type discussed above was that of Infeld and Wallace (IW). jl^] 
This derivation applied the EIH procedure to the coupled Einstein-Maxwell equations and was thus able to avoid the 
difhculties of the other derivations. Strangly, this paper is almost never refered to in the literature. However, in spite 
of its superiority over other derivations, the IW derivation did not overcome all of the difficulties of the ALD equation. 
IW gave seperate derivations using half-advanced plus half-retarded and pure retarded potentials and argued that the 
former was more "natural" even though there was no radiation damping in this case! In their derivations IW focused 
on one charge, which they took to be momentarily at rest, in an external field. After obtaining their conditions of 
motion for it from EIH type surface integrals they then performed a Lorentz boost to obtain conditions for arbitrary 
velocities. Unfortunately, dynamical effects involving accelerations can be missed in this way and indeed, such effects 
are not present in the IW conditions. They also assumed, as did EIH, that they had derived equations of motion 
rather than asymptotic conditions and hence were unable to rule out run-away solutions. And finally, theirs was a 
formal expansion in powers of a parameter A which was in the end set equal to unity. 

In what follows I will present the details of the derivation of the conditions imposed on the motion of the sources 
of electromagnetic and gravitational fields by applying the EIH method to the coupled Einstein-Maxwell equations 
of general relativity. In the next section I will discuss briefly earlied derivations of the radiation reaction force with 
emphasis on the Dirac derivation which is nearest in spirit to the EIH method so that one can see the essential 
differences in the two approaches. Section III will present a derivation of the EIH surface integrals and Section IV 
will outline the approximation procedures needed to evaluate these integrals. Section V contains the details of the 
derivation and the paper concludes with a brief summary. 

II. CLASSICAL DERIVATIONS AND THEIR PROBLEMS 

By classical I mean here derivations of the electromagnetic radiation reaction force that do not use the general theory 
of relativity. Most of these derivations attempt to calculate in one way or another a "self" -force of a charged particle 
acting on itself. The first derivations of this kind were given in 1903 by Abraham and 1904 by Lorentz. |jl^ These 
early derivations made use of a charge model consisting of a rigid, spherically symmetric charge distribution p(x) of 
radius a and calculated this force by summing up the interactions of pairs of the infinitesimal elements comprising the 
charge. Because of the assumption of rigidity the results of these calculations were only applicable to non-relativistic 
motions. The result of these calculations was a series expansion of the form 
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F,elf=^F„ (2.1) 



where 



and 



4 U 

Fo = -3^v, (2.2) 
F.^^v, (2.3) 

Here v is the velocity of the charge, e its total charge and U its so-called self-energy which is given by 

It is this term which gives most of the trouble in these derivations since it is proportional to e^/c^a. In order to 
eliminate the higher order terms in the expression for Fgcif one would like to take the limit a — > cxd. But then of course 
U diverges. There is however, in my view, an even more serious problem with these derivations. In order to evaluate 
Fscif one must make use of the Lorentz force law. But that use does not in and of itself justify inserting it into the 
right hand side of the Newtonian equations of motion, a procedure I would call proof by naming. Nor can one be 
sure that the form of Newton's inertial force, mo v or its special relativistic generalization, is valid when dealing with 
charged particles. If one does proceed in this manner, then one is forced to take for tuq, the expression 

'ma^m-- — , (2.6) 
o c 

where m is the physical (finite) mass of the charge, in order to obtain a finite result in the limit a — > cx). In this limit 
the bare mass toq is of course negatively infinite. 

In order to avoid the non-relativistic restrictions of the AL model, Feynman Q suggested using a cutoff procedure 
in which the Fourier representation of the retarded Green function is altered by making the substitution 

1 _ 1 (2.7) 

where A is the cutoff factor and is ultimately allowed to become infinite. Using this modified Green function, one 
can calculate the self-force of a point charge on itself by evaluating the Lorentz force at the location of the charge. 
1^ The result is finite and is the relativistic extension of the AL self force with the charge radius a replaced by 1/A.. 
Aside from this feature, this derivation suffers from the same defects of the AL derivation including the infinite mass 
rcnormalization. 

A different approach was taken by Dirac |0 that has some of the features of the EIH derivation and is reviewed 
here so that the reader can compare the two derivations. Dirac considered point charges and defined what he called 
the electromagnetic radiation field F^J^^ to be 

^rad = ^^^^ - F^'"' (2-8) 

where and F~^^'^ are respectively the retarded and advanced fields of a such a charge. While these latter fields 

are infinite on the charges trajectory, the radiation field is not and is given there by 

T-nUl/ 4:6 ,...fj. . I, ...iy.fi-. /n n\ 

^rad = y ^ - ^ ^ ) ■ (2-9) 

where z'^{t) are the coordinates of the charge expressed as a function of its proper time t and where dots over a 
quantity indicate differentiation with respect to r. Dirac argued that his definition of the radiation field was physically 
reasonable since, at large distances from the charge and at the correspondingly large times after an acceleration takes 
place, the advanced field will be zero. He then assumed the existence of an energy-stress tensor 
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matter 



rpflV 



(2.10) 



where T^attor ^ matter energy-stress tensor localized on the particles trajectory and 



(2.11) 



is the electromagnetic energy-stress tensor. In these equations rjf^i, = diag(l, -1,-1,-1) is the Minkowski metric tensor 
and T]^'^ is its inverse , Greek indicies take the values 0,...3, and the Einstein summation convention is used throughout. 
The total energy-stress tensor is assumed to be conserved so that 



0, 



(2.12) 



where the comma notation is used to denote differentiation, that is ^ = d^t^- 

Dirac proceeds by integrating this conservation law over the four-volume contained within a thin tube surrounding 
a piece of the charge's trajectory and converting the result into an integral over the surface of the tube by means of 
Gauss's theorem to obtain 



(2.13) 



where ai and (72 refer to the top and bottom caps respectively of the tube and S refers to its wall. Since, by assumption, 
the support of the matter energy-stress tensor is the charge trajectory it does not contribute to the integral along the 
wall. On it Dirac evaluates using the total field F/^^ given by 



^tot 



(2.14) 



where F-^^^ is the field incident on the charge. The wall is such that its intersection with each space-like plane whose 
normal is parallel to at the point of its intersection with the charge trajectory is a small sphere of radius e. In 
order that this construction leads to a well-defined wall it is of course necessary that e be less than the curvature 
of the trajectory which is inversely proportional to the acceleration of the charge. Finally the field on each of these 
spheres is evaluated in terms of the coordinates of its instantaneous center. This evaluation cannot be given exactly 
since F'^'^" depends on the retarded position of the charge which is given by an expansion in e about its instantaneous 
position. The result of these calculations is that 
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dT 



(2.15) 



where Q and Q' correspond to the points of the trajectory that intersect the top and bottom caps respectively and 
where 



r 



Ft.- 



F- 



(2.16) 



11) above for 



The 0{e) term in (|1|) is in fact an infinite series in e similar in form to the Fn in the expression 
fscif which Dirac does not evaluate. 

In order to complete the derivation it is of course necessary to evaluate the cap contributions in ( |2.13| ). But here 
one encounters two serious difficulties: one does not know T^attor '^^ initio and worse, the contribution to the caps 
from is divergent. As a consequence Dirac was forced to assume the form of the contributions from these caps. 
Since the final answer in any case must be independent of e these contributions must cancel the e dependent terms in 
the integral in (2.15) and not, as Dirac asserts, so that the integrals have a definite form in the limit e — > 0. There is, 



in fact, no need to take this limit nor would it make any mathematical sense to do so. These contributions must also 
be consistent with the requirement, which follows from (2.13), that the line integral in ( 2.15D depend only on the end 
points Q and Q' and hence the integrand in this integral be a perfect differential. Neglecting the e dependent terms 
which in all cases must cancel with contributions from the caps this latter condition tells us that 



/.^ =i? 



(2.17) 



where is undetermined except for the requirement that i^^B = which follows from the antisymmetry of f^" in 
fj, and ly. The simplest assumption that satisfies this requirement is that = m z^although, as Dirac pointed out, 
there are other possibilities. With this assumptionn the equations of motion become 
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m 2 



[z^ + rz^z''] + e z, (2.18) 
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which Dirac took to be exact. However, buried within the assumption for is some kind of infinite mass renormal- 
ization to compensate for the divergent self-energy contribution of T^'^ to the cap integration. If, for example, one 
took, for the matter energy-stress tensor the Minkowski tensor 

TZuer = rno f 5\x - z(r)) z^ z" dr (2.19) 



where 5 {x) is the four-dimensional Dirac 5-function, the bare mass mp would contain a negative infinite component 



similar to the ttiq in (2.6). As the reader will see in what follows, the EIH surface integral method circumvents this 
problem by in effect replacing the volume integrals in the cap contributions by surface integrals which are all finite. 
Finally, this equation, like its counterparts obtained from the use of one or the other radiation reaction force, possesses 
run-away solutions, that is, solutions for which the accelerations increases exponentially with time. 

III. EIH SURFACE INTEGRALS 

The derivation of the conditions of motion makes use of the fact that the Einstein-Maxwell field equations for 
the gravitational field and the electromagnetic field F^^'^ contain three dimensional curls. Therefore, when these 
equations are integrated over any closed 2-surface these curl terms vanish. When in particular such a surface encloses 
a source of these fields the requirement that the remaining surface terms must be surface independent and vanish 
leads to conditions on the motion of this source. To facilitate the construction of these surface integrals it is useful to 
write the Einstein-Maxwell field equations in the form (I use units in which G — c — 1^ Latin indices run from 1 to 3 
and Greek indices from to 3.) 

U^'^'P^p = 6'^'' (3.1) 

where 

^^.p ^ _^ppu ^ J_{(_5)(5M^'^P- _ g^Pg-)}^ , (3.2) 

Q^"" = {-9){TZ + tll). (3.3) 

(V^F^'^),, = 0, (3.4) 

and 

F[tJLV,p\ = Q. (3.5) 

Her e g = det(ffmy), t^^^ is the Landau-Lifshitz pseudotensor T^i^ is the electromagnetic energy-stress tensor given 
by ( ^.ll| ) (in what follows I will omit the subscript em to simplify the notation), indices are raised and lowered using 
g^i/ in the usual manner and brackets around a set of indices indicates that the indices are to be antisymmetrized. 
Because of the antisymmetry of F^'^ in ji and v and JJ^^'^f in v and p, it follows that and U^J''' are three 



dimensional curls. As a consequence, integration of (3.1) over a closed spatial 2-surface in a t =const. hypersurface 
gives the result 

{Uf^^-^^o-Q^nrdS ^0 (3.6) 
where is an outward pointing unit surface normal. Similarly, we get from equation (2.4) 

{V^F''°)^onrdS = 0. (3.7) 



It is these last two equations that are used to obtain conditions of motion for our charged system. They come 
from noting that, since these equations must hold on any two-surface, the sums of the contributions to the integrals 
that are surface independent must by themselves be zero. (The surface dependent terms must in all cases vanish 
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identically.) In practice one chooses the surfaces to be spheres centered on the sources and sets to zero those terms 
that do not vanish identically. (If there are no sources, so that the fields are divergence free inside the surfaces, then 
all the integrals vanish directly as a consequence of the field equations.) 

Because one employs in the derivation only surface integrals on which the fields are every where finite, it follows 
that there are no infinite integrals that must be renormalized away. If we compare (3.6) with (2.13) we see that the 
divergent cap terms in the latter equation get replaced by the finite surface integrals of C/^''°^o and {—g)t^]^. It is the 
first of these integrals in fact that contribute the inertial force terms to the conditions of motion, obviating the need 
to postualte them. Furthermore, since the source region is contained entirely within the 2-surfaces of integration in 
(3.6) and (3.7), one need only know the fields on these surfaces and not inside. All that is required is that each surface 
contain only one source, that their radii are sufficiently large that the fields on them are weak and small compared 
to the distance between sources. Thus one can equally well deal with compact sources with strong internal gravity, 
e.g., neutron stars and even black holes as with sources whose fields are weak everywhere. In particular, the problems 
associated with the use of point sources discussed in Section II do no occur in the EIH procedure. 



IV. APPROXIMATION PROCEDURES 



The evaluation of the surface integrals in (3.6) and (3.7) requires a knowledge of the fields on these surfaces and 



this can only be gotten by employing some kind of perturbation scheme. The use of any such scheme requires that 
there be one or more small parameters associated with the system. In the present work one of these is a dimensionless 
'slowness' parameter e that is of the order of the ratio of the light travel time across the system of sources to a time 
scale associated with their motion, e.g. an orbital period. It is also of the order of the ratio of the size of the system 
to the wavelength of the radiation emitted by it and also the ratio of the velocities of the sources to the velocity of 
light. Its precise value is fixed by the parameters that characterize the source system and its motion. It should be 
emphasized that this parameter is not an arbitrary small parameter as EIH took it to be nor can it be taken to be 
the inverse of the velocity of light as many authors have done. In addition to assuming that e is small compared to 
unity, I will assume that the relative strengths of the gravitational and electromagnetic interactions 6 =mass/charge 
is small compared to unity and also to e for all charges in the system. In this work I have only kept terms of order 
5^ and I have also scaled masses and charges so that ttia = 0{5^€^) and qA — 0{5e^). These restrictions have the 
effect that the operation of raising and lowering indices on all tensorial quantities is accomplished with the help of 
the flat-space Minkowski tensor. 



A. Matched Asymptotic Expansions 

The perturbation methods employed here are the singular perturbation techniques used in AI to derive approximate 
conditions of motion for gravitationally interacting uncharged particles. The first of these is the method of matched 
asymptotic expansions which was first used by Laplace in 1805 (see Ref. [|l8| for this and other references to this 
method.) In general this method seeks to develop approximate solutions to a system of equations in different regions 
of the range of the independent variables using expansions appropriate to each region The two expansions are then 
'matched' in an overlap region where both expansions are valid. In order to obtain approximate solutions to the 
Einstein-Maxwell equations it is in principle necessary to use four overlapping regions in all corresponding to two 
different matchings, one in time and one in space. 

The time matching treats the region in the vicinity of an initial value surface as a boundary layer called here the 
fast zone where the solution depends on the initial value data. The solution in this region is matched to an asymptotic 
region in the future of the initial data surface t = const, called the slow zone. The undetermined parameters and 
functions appearing in the solution in this latter region are then determined in terms of the initial value data by 
matching to the fast zone. The slow zone can be thought of as the region in which all initial transients generated by 
the initial data have died off. 

The other matching is the one used by Burke. [po| It matches solutions in two overlapping spatial regions consisting 
of an inner or induction zone and an outer or radiation zone. What is important to us is that, in the lowest order of 
approximation the solution in the outer, slow zone is a pure outgoing solution, that is, a function of a null coordinate u 
and the spatial coordinates r. (In higher orders of approximation this region will contain both outgoing and incoming 
radiation due to backscattering.) It should be emphasized that this form of the solution is not imposed arbitrarily as 
it usually is in other treatments. Rather, it is a consequence of the fact that solutions of the wave equation have the 
property that so long as the energy of the field on the initial-value surface is finite then in the asymptotic future the 
field is purely outgoing. pl| 
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As we shall see, it is the matching to this outgoing solution that is responsible for the phenomena of radiation 
damping and the appearance of irreversible terms in the conditions of motion to be derived here. Thus we see that 
irreversibility in the conditions of motion is a consequence of the finite energy condition imposed on the initial fields 
and not, as IW contended, because of an arbitrary choice in the use of retarded or advanced potentials. One could 
choose initial conditions on the fields that would lead to anti-damping for all times but this would violate the finite 
energy condition mentioned above. While the final result of the two approaches appears to be the same, the one used 
here is in accord with those used in other areas of physics where irreversibility arises in the asymptotic future as a 
consequence of restrictions on initial conditions. 

Although in principle one could construct solutions in the fast-time zone to obtain conditions of motion in this 
zone, I have not explicitly done so in this work. In general, the motion of the sources in the fast-time zone will depend 
in a very complicated way on the initial values not only of the source coordinates but those of the fields as well and 
will almost certainly not satisfy any unique set of equations of motion. Without this matching the conditions on the 
source motion obtained by the other matchings will then only be valid in the asymptotic future and are therefore not 
equations of motion in the usually accepted Newtonian sense. In particular, they cannot be used to solve initial value 
problems. I will defer a discussion of the implications of this state of affairs until after I outline a derivation of the 
conditions of motion accurate to radiation damping order. 



B. Multiple Time Scales 



In order to construct approximate expansions in the two spatial domains in the slow-motion domain it is necessary 
to make use of the method of multiple time scales. In their original work EIH employed only a rudimentary multiple 
time scale formalism. As a consequence they were forced to introduced fictitious dipole terms at each order of their 
expansion in A and in the end require that the sum of all these dipoles vanish. This procedure appeared somewhat 
arbitrary and may have cast some doubt on the validity of their results. In fact, when the multiple time scale formalism 
is used no such dipole terms are needed. 

As the name implies, the multiple time scale formalism is used when one has to deal with systems in which 
phenomena occur on different time scales, e.g., oscillation and damping. Not only is it useful in separating out these 
different types of motion but it is used extensively to avoid the appearance of secular growth and nonuniformities in 
various approximation schemes. In its application one assumes that the dependent variables of a system of equations 
depend on the time variable t via a sequence of times, to,ti, ... where usually, but not always, to = t and tn — e"i, 
e being a small dimensionless parameter in the problem. In this case, which is the one that applies in this work, a 
derivative with respect to t is replaced by 

dt ^5]e"5t„. (4.1) 

The dependence on these times is then fixed by the requirement that the terms responsible for secular growth vanish. 
(See AI and Ref. iQ for simple examples of the application of this method.) 

The time scales required in general relativity will depend upon the problem being considered. If we are dealing 
with uncharged particles and if t is the fastest time, characterized by the light crossing time of the system, then 
et will characterize the Newtonian time scale, e^t the post-Newtonian effects such as perihelion advance, e^t the 
post-post-Newtonian effects etc.. Gravitational radiation is characterized by e^t, sometimes referred to as the 2^- 
post-Newtonian time scale (see Reference 11 for details.) When one is dealing with charged sources there are at least 
two different small parameters that characterize the system, e and S = m/q, and hence the possibility of additional 
time scales arises. However. I have assumed that 5 <C e ^ 1 so as to avoid having to deal with both electromagnetic 
and gravitational effects simultaneously. Consequently, ti = et will characterize the Coulombian interaction time 
scale, i3 — t^t the Amperian or post-Coulombian time scale and radiation effects will be characterized by a ^4 = e^t 
time scale in lowest order. [Esi 



V. CONDITIONS OF MOTION 



The constr uct ion of app roximate solutions of the field equations (3.1) and (3.4) and their use in evaluating the 



integrals in (B.6) and (3.7) follows the pattern developed in AI. (Here I will use the vector 4-potential where 
F^n — A^,^ — A^^i,. As in AI, I assume the existence of a coordinate map such that g = 33'^" can be expanded up 
to the order of accuracy needed here in an asymptotic series of the form 
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(5.1) 



where I have assumed that the /i^J:^ depend on the multiple times through their dependence on the source coordi- 
nates x^(ti, ^3, . . . , e, (5), charges qA{ti, ^3, ■ ■ • , e, 5) and masses m^(ii, is, . . . , e, 5) with A labeling the sources. One 
should keep in mind that q^'^ cannot be expanded indefinitely in a simple power series since it is not analytic in e, 



reflecting the fact that one encounters e" ln(e) terms in higher orders of the approximation 
do not appear at the level of accuracy used in this work. In a like manner I assume that t 
expanded as 



24[ | However, these terms 
re 4-potential can be 



E 



e"(5"MJ^,„(r,XA, qA, ruA) 



(5.2) 



(Recall that, because of the smallness of 5, indices are raised and lowered here using the Minkowski tensor so that, 
e.g., A^" = Ty'^M,.) 

starts with 0{5^) terms while the 



Because of th e scalings chosen for our sources it follows that the expansion (5.1 



expansion (5^) starts with 0{5) terms. Since one only needs terms of this order to obtain an 0{5'^) overall accuracy 
in the conditions of motion I will drop the index m in the above expansions in what follows. Also, in this order the 
dependence of the source variables x^i , qA, mA on S can be ignored so that from this point on I will ignore entirely 
all 6 dependencies. Then, again becau se o f the scaling employed, the lowest order in e fields are A2 and /i™. 

In addition to the field expansions (5J) and ( ^.2| ), it will also necessary to make certain assumptions concerning 
the dependence of the source variables qA, mA, and xa on e. To the level of accuracy needed here I will assume that 
the latter two variables can be expanded as an asymptotic series in powers of e of the form 



TUA 



e*mA4 



(5.3) 



and 



X^ XAO 



e^XA2 



(5.4) 



where the dots indicate terms that are at least an order of magnitude in e smaller that the terms retained. It should 
be emphasized that the powers of e used in this expansion is not arbitrary but dictated by approximation itself. As 
we shall see, the masses ttia will be automatically expanded as we proceed in the approximation while the charges qA 
are absolute constants of the motion and that there is no need to e xpand them at all. 

To obtain equations for the fields /i^'' and from (3.1) and (3.4) I will impose the deDonder coordinate conditions 



and the Lorentz gauge condition 







0. 



(5.5) 



(5.6) 



Here too problems arise with the use of these conditions in higher orders of the approximation and they must be 
modified to eliminate the appearance ln(r) terms in these higher orders. p4[ 

At this point a difficulty arises. Usually when one substitutes an expansion of the dependent variables in powers 
of a small parameter into a system of equations one simply equates to zero the coefficients of the var iou s powe rs o f 
this small parameter. However, in the present case the coefficients /i^^ and Al^ in the expansions (5T) and ( |5.2| ) 
themselves depend on e through the dependence of the source variables on e. Therefore one cannot simply equate to 
zero the coefficients of the powers of e when these expansions are substituted into the field equations. This difficulty 
seems to have been overlooked in all previous derivations of equations of motion in general relati vity includ ing AI as 
far as I know. It can however be dealt with in the foll owin g manner: Substitute the expansions ( |5.l| ) and (5.2) into 
the field equations and set e = 0. It then follows from (3T) that in the source- free weak- field region/i2*'|c=o satisfies 



2, OOi 



0. 



(5.7) 



(5^ 



However, since /i^*^ depends on e solely through its dependence on the source variables it follows that it must satisfy 
for all values of e ^ 1 . Since at each stage of the approximation the unknown field coefficients in the expansions 



and (5.2) enter the field equations linearly it follows that this argument holds quite generally and so allows us 



to proceed to find these coefficients without regard to the dependence of the source positions on e. 
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The requirement that the sources are characterized only by their charge and mass is equivalent to requiring that, 
in the absence of other sources, the field of each source must be spherically symmetric 
is given by 



It follows therefore that h^*^ 



(5.8) 



where the sum is over all sources in the system, i.e., over A = 1 to N ,r\ = r"^ — x\ where x\ is the «th coordinate 
of the Ath source and where, for convenience, we will take the center of mass of our system to be at the origin of 
coordinates so that ^ mAx\ — 0. (In the case of a charged black hole, x\ is the ef fect ive coordinate of its center as 
determined by its exterior vacuum solution.) In a like manner the field equations (3^) allow us to conclude that in 
the source-free region A2 satisfies 



2 



and so is given, for spherically symmetric sources, by 







Al=Y,(lAlrA. 



(5.9) 



(5.10) 



It cannot be emphasized too strongly that these solutions hold only in the exterior of the field sources and are not 
in any sense 'singular'. There are in fact no singular fields in the usual sense in general relativity as there are in special 
relativistic theories for example. The sources can be either compact matter sources or black holes which are modeled 
by our ass umpti on th at they have no internal dynamics and are characterized solely by their mass and charge. The 
fields (5^) and ( 5.10| ) are exterior fields and by Birkhoff's theorem are unique. 

It should also be explained why the coordinate ta appears in the above solutions and not it's zeroth approximation. 
Indeed, if we were to substitute the expansion (5.4) into these solutions we could then expand them about e = 0, 
retain in them only the lowest order terms and add the remaining higher order terms to the higher order terms in 
the expansions (5.1) and (5.5). In effect, by proceeding in this manner we have effectively summed partial series that 
would have resulted had we not done so. By keeping all of these terms in the lowest order solutions (and we will 
continue this practice with all of the higher order terms) we accomplish two things: Firstly, it greatly simplifies the 
calculation due to the effe ctiv e sum min g of numerous partial series. But secondly, and most importantly it enables us 
to center the surfaces in (3.6) and (3.7) on the true centers of the sources rather than on their approximate centers. 
Although of course one can use any surfaces one likes in these equations it again greatly simplifies the resulting 
calculations to use the true source centers since doing otherwise results in spurious multipoles appearing in the higher 
order approximations. 



A. Coulombian conditions of motion 



With the help of the fields and A^ found above we are now able to determine the time dependence of the 
charges an d m asses appearing in t he fi elds given above. To do so one needs to evaluate the surface integrals that 
appear in (3^) with ^ = and in (3/7) respectively using these approximate solutions. For convenience one chooses 
spherical surfaces centered on the sources with radii large compared to their physical or Schwarzschild radii in the 
case of black holes but small compared to the distance between them. These restrictions on the size of these surfaces 
are necessary to insure that the weak field approximation employed here is valid in their neighborhoods. In carrying 
out these evaluations one finds two types of terms - those that depend on the radii of the surfaces and those that 
don't. Since the integrals as a whole must be independent of these radii it follows that the overall coefficients of the 
various powers of the radii must vanish identically. It also follows that the radii independent terms must also vanish 
since the integral as a whole must vanish. However, these terms will not vanish identically. The requirement that 
they do vanish imposes conditions on the motion of the sources as EIH showed. As a consequence, one need only 
pick out those terms in the integrands that are proportional to the inverse squares of the surface radii and direction 
independent and set their coeffici ents equal to zero to obtain these conditions. 

Substitute then the expansion (5T) into (3^), divide by and set e equal to zero. In this way one finds that 



dtimA2 = 0. 



(5.11) 



The surface integral in equation (3/7) can be evaluated using A2 given by ( 5.10| ) without the need to set e equal 
to zero since it is linear in A^ and since higher order terms in the expansion ( ^.2[ ) will not contribute to the integral 
because of their angular dependence. As a consequence we obtain the result that 
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dtqA = 0. 



(5.12) 



The qA are thus absolute constants of the motion. Note that a sim ilar argument does not pertain in the case of the 
masses because, in higher orders of the approximation, in (3.6) will in general be non-zero. 

The next quantity needed to obtain the lowest order conditions of motion is the lowest order contribution to the 



field h . In order to satisfy the deDonder conditions 5.5 we see that these contributions must be of order e'^ since 
is of this order and differentiation with respect to a spatial coordinate does not change the order of a quantity. 
It then fo llow s from the field equations and the deDonder conditions, using reasoning similar to that which led to the 
solution (^), 



that ft.g'' satisfies 







(5.13) 



and therefore is given by 



mA2 XAr/rA 



(5.14) 



where a dot over a quantity denotes differentiation with respect to ti . 

Equation ( |3.6[ ) with /i = s can now be used in conjunction with the above fields to obtain the lowest, Coulombian 
conditions of motion. For this purpose we need to know T^"^ to 0{€^) and JJ^"^^ to 0{e^). The integral over the spherical 
surface surrounding the Aih source appearing in ( |3.6| ) will vanish provided the overall coefficient of the terms in the 
integrand proportional to 1/r^ that remain after dividing by e"* and setting e equal to zero vanishes. This will be the 
case if 



"^A2 i^o^ IaF^a 



(5.15) 



where F^^ is the lowest order electric field at the location of the Aih. source due to all the other sources and is given 
by 



^2A 



E 



^AB 



"ABO 



(5.16) 



where = — and the prime on the sum indicates that it is over all B ^ A. Note that, at this order of 
approximation, our conditions of motion involve only the lowest order contributions to the masses and positions of 
the sources. However, unlike in the case of the fields, we cannot conclude that these equations are valid for non-zero 
values of e, i.e., we canno t drop the subscripts in ( 5.15| ) and (5.16). It must be emphasized that the excl usio n of 
A from the sum in (5.16) is not arbitrary but rather is a direct consequence of evaluating the integrals in (3.6). In 
fact it could not be otherwise since all integrals are over surfaces and involve only finite quantities there. It must also 
be emphasized that the inertial terms in these equa tions are not inserted by hand by appealing to Newton's laws of 
motion but are already implicitly contained in (3.6). 



B. Post-Coulombian conditions of motion 



To o btain the next, post-Coulombian, conditions on the motion of the sources one needs to evaluate the integrals 
in ( |3.6|) with /i = s to an accuracy of O(e^). The Maxwell stress-energy tensor T'"" appearing i n th e integrand can 
be evaluated to this order from a knowledge of and Al. It follows from the field equations (3.4) and the gauge 
condition ( p^ ) that A^ satisfies 

V'^Al = 0. (5.17) 
Therefore, as a consequence of the gauge condition, it must have the form 

^3 = E'^'4 Ita- (5.18) 
These same equations and gauge condition also lead to the equation 

^^Al- A^= 0. (5.19) 

It follows that Ajj is given by 
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QA rA 



z ta 



+ Xa n. 



Xa\. 



(5.20) 



where n\ = r\/rA- No homogeneous solution of equation ( 5.17 ) is included here since, by assumption, it would have 
to be of the form (5.10) and would therefore have the effect of merely redefining the charges qA appearing in these 
solutions. This expression for A^^ appears to contain a term involving the acceleration of the ATh. source. However 
the conditions (5.15) can be used to replace the acceleration factors with Coulomb interaction terms without affecting 

the overall accuracy of the expression, i.e., it will still be accurate to 0{€^). 

With these results in hand the evaluation of the second surface integral in ( p.6[ ) over a spherical surface surrounding 
the ATh source to 0{e^) is straightforward. Again, one only needs to determine that part of the integral that is 
independent of rA- One finds in this way that 



dS^ = XAr 



^3A 



„ TpmO 
lAi'iA 



(5.21) 



In this expression F^^^ and F4 ™ are again the fields at the location of the Ath source produced by all of the other 
sources. The first of these fields is given by 



FsArs = ^ {A3As,r ~ ^3Ar,s) QB -^i^B ^AB' ^B ^AB )• 



(5.22) 



''AB 



Likewise the second field is given by 



FiAOr — (AsArfi ~ ^4A0,r) 



2 



{ r f ' ^ s \1 
— Xj^g{Xg Xab) + 

^AB 



"■AB 



-{x 



AB^AB ^B ~^AB 



■ 2 

X R 



XAB 



(5.23) 



where again the acceleration factors appearing here are to be replaced by Coulombian interactions with the help of 

(Hi- 

It appears that somethi ng is wrong with our calculation, since, among other things we only get 2/3 of the expected 
Amperian interaction in (5.21) and the sign of the second term on the right hand side of t his equations appears to 
be wrong. We are led to this conclusion however only if we try to interpret equation (|3.6|) itself as an equation of 
motion so that we might be led to expect that <f U^"^ 



surface over which the integral is performed and § 0'" 



. dS is the 4-momentum of the source contained within the 
r dS is the 4- force acting on it. In fact they are not. What 
happens is that part of the first integral contributes the missing 1/3 of the Amperian interaction even though the 
electromagnetic field does not appear explicitl y in its definition. For its evaluation we need to know {/^"O .y^,l^ich, 
when use is made of the deDonder conditions (5.7), is given by 



7-7-mnO 

^^5 — 



(5.24) 



It follows from the field equations (3.1) that /i™" satisfies 



2j mO 



..mO 



-AF^^F™ = 



(5.25) 



while ft™" satisfies 



•Or 
2 



4 pOm pOn ^ Q_ 



(5.26) 



Thus we see that the electromagnetic field acts, at this order, as a source of the gravitational field and it is these 
contributions that in turn supply the missing 1/3 of the Amperian force. 
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It is not possible to obtain closed form solutions for either of the above equations because of the terms that are 



quadratic in the fields. However, all that is needed is to determine those parts of the unknown fields /i™" and /i™" 
in the neighborhood of a particular source that will yield surface independent terms in order to evaluate the surface 
integrals surrounding this source. The needed terms are found following procedures developed by EIH and in AI 



and are given in the appendix. To obtain the post-Coulombian conditions evaluate the surface integrals in (3.6) with 
fi — s, take the second derivative with respect to e and set e equal to zero. The result of these operations yields the 
conditions 

mA2{Xj^2 +^ ^A0.t3 + ^AO^AO + ^AO^AQ^AOf = 



QA{^x';^ggF^QX^J^^2 



J^4A0 



1 .2 
- X 

2 



AO ^2A 



'I 11,1 L ™ I 

SAO ■''AOS 



(5.27) 



As in all multiple time approximations, the dependence of a;™Q on is in this equation is determined by the requirement 
that the totality of the terms that would contribu te to the secular growth of niust vanish. 

These conditions, together with the conditions (5.15), are identical in form to the first two slow motion approxima- 
tions one would obtain from the special relativistic equations of motion for a charged particle of mass mA2 and charge 
qa moving in the field of a collection of other charges using a multiple time expansion. The second and third terms 
on the left side of this equation as well as the second term on the right side are equivalent to terms in the special 
relativistic equations that arise from replacing derivatives with respect to proper time by derivatives with respect to 
coordinate time and expanding the factors of 1/^1 — (e i )2 that enter as a consequence of the replacement. The 
first term on the right side is equivalent to a term that is due to the effect of retardation while the last term is 
the Amperian force. It should be emphasized again however that the conditions derived here differ from the special 
relativistic equations in a fundamental way in that they are only asymptotic conditions and are not required to be 
valid at all times. They are also derived from a more fundamental theory without having to be postulated ab initio. 

Th e reader should notice that, without the multiple time formalism, the first two terms on th e left side of equation 
( 5.27 ) would be absent and we would arrive at two totally different equations, ( 5.15| ) and ( 5.27 ), for the accelerations 



and hence a contradiction. In an attempt to avoid this inconsistency some earlier treatments expanded the 
coordinates xa in a series in e. However, this procedure by itself will fail without the multiple time formalism since 
secular growth, as in the case of periastro n adv ance, will ultimately lead to a nonuniform expansion. 

In arriving at our conditions of motion (5.27) we have consistently ignored surface dependent contributions to our 
surface integrals since the sum of these contributions must v anish identically. One such contribution does deserve 
special comment however. In evaluating the surface integral (5.21) one finds a contribution | /t a- As ide from a 



difference in the numerical coefficient, this is the same contribution appearing in the line intcrgral in ( 2.15 ) f oun d by 
Dirac. Here it is exactly canceled by a corresponding term coming from first term in the surface integral in (3.6). 



C. Radiation reaction 

The effects of radiation reaction can be derived using the matched asymptotic expansion method of Burke pO[ | and 
its extension to the surface integral method in AI coupled with the method of multiple time scales. Since Burke did 
not use the surface integral method he had to make an assumption equivalent to assuming the validity of the geodesic 
equations in order to compute the radiation reaction force. Further, the fact that he did not use multiple time scales 
meant that his result cannot be considered to be part of a consistent approximation scheme even though he did obtain 
the "standard" form of the gravitational reaction force. 

To determine the electromagnetic reaction force I will proceed as in AI by constructing an outer, slow zone solution 
and match it to an appropriate inner solution in the neighborhood of a given source. To this end, following Burke, 
I introduce outer coordinates x'^ = (ei, er). The fields A'^{x) = A^^{x) are again expanded in what is hopefully 
an asymptotic series in e where now the coefficients in the expansion are functions of the outer coordinates. If one 
imposes the Lorentz gauge condition it then follows from the field equations that, to order S"^, A'^ satisfies the fiat-space 
homogeneous wave equation 

(V^ - 8^)1" = 0. (5.28) 

Since the inner expansion of the outer solutions of ( [3.2§| ) must match to the outer expansions of the inner solutions 
we have already constructed we need these latter expansions. They are obtained by letting r ^ oo while holding e 
fixed. Keeping only the lowest order nontrivial terms needed for this purpose in this limit yields the results 
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(5.29) 



where 



qA, 



(5.30) 



and 



P_ 

r 



(5.31) 



We see from these outer expansions that the lowest-order contributions to A° must contain terms with zero- and 
first-order spherical harmonics as factors while those to must contain zero-order harmonics. 

While the spherical harmonic structure of the outer fields restricts their form it does not fix it - any linear com- 
bination of outgoing and incoming wave with the appropriate harmonic dependence will do. In most treatments of 
radiation it is assumed that only the outgoing wave is present in the wave zone and indeed such was assumed both 
by Burke and by me in our respective previous works. However, as I later showed, p]| if the energy contained in 
a wave field on some initial value surface is finite, then in the asymptotic future of this surface the field will be a 
pure retarded wave. (If one is willing to violate this condition then it is possible to have incoming waves for as long 
as one wishes into the future.) One could also have solved the wave equation as a final value problem. If one now 
required that the wave energy on a final value surface be finite then the solution would be a pure advanced wave in 
the asymptotic past. As a consequence, we see that the choice here of the retarded wave solution to (5.28) is not 
arbitrary. At the same time it means that any results that depend on this choice such as the about to be derived 
radiation reaction force will only by valid in the asymptotic future. 

Retarded waves are functions oi u = t — r and r. Consequently the ones that have the right angular dependence 
are 



and 



A'' = 



9r{u, e) 



(5.32) 



(5.33) 



where the prime in equation (5.32) denotes differentiation with respect to u. The as yet arbitrary functions /r(u) and 
gr{u) appearing in these solutions are not independent. They m ust b e so chosen that the field s A^ and A"^ satisfy 
the Lorentz gauge condition used in obtaining the wave equation ( ^.2j ) from the field equations (3^). The latter will 
be the case provided 



(5.34) 



To determine these arbitrary functions completely we must match the inner expansions of the outer solutions ( 5.32| ) 
and (5.33) to the outer expansions of the inner solutions ( 3.2g| ) and (5.31). These outer expansions are obtained by 
letting r — !■ while holding e fixed and are given by 



h 



frjet) 
(er)2 



I 'fr i^t) + ^ 'f\. iet)er + 



and 



A'' 



e r 



(5.35) 



(5.36) 



where again a dot over a quantity denotes differentiation with respect to ii = e t. The inner and outer expansion will 
match if we take 



b^e^Q and f riet) = e^P'' (et). 



(5.37) 



In the inner expansion ( ^.35| ) the first two terms match to the first two terms in the outer expansion of A n, the 
third term matches to the first term in the outer expansion of A^ and the first term of the inner expansion ( 5.36| ) 
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matches to th e first term in the outer e xpan sion of A^. However, the time-odd terms in these expansions, the fourth 
term in ( 3.35| ) and the second term in ( 5.36| ), do not match to terms in either of these outer expansions. Since both 
of these terms are sohitions of the homogeneous Laplace equation a matching can be achieved by adding solutions 
of the Laplace equation to our inner solutions whose outer expansions match to these time-odd terms. The required 
terms are 



1 



(5.38) 



and 



(5.39) 



The above potentials can now be used to evaluate additional surface independent contributions to ( |3.q ) with fj, = s. 
These lowest order contributions to the second integral in this equation are C(e^) and, for a sphere surrounding the 
Ath source, are given by 



T'''n\ dSA 



E 



(5.40) 



Note in particular that in this expression the sum is over all of the sources including the Ath source itself and thus 
includes, for the first time in these calculations, the effect of a self action. We shall see that this self action is similar 
to the radiation reaction in the ALD equations. 



by 



There is also an O(e^) surface independent contribution to the first integral in equation (3.6) and is given simply 



W^^rfj, dSA 



2771 AO dt. 



(5.41) 



where — e^t. The 0{e^) conditions of motion arc then gotten by dividing equation (3.6) by e^, differentiating five 
times with respect to e and finally setting e equal to zero. The result of these operations is 



2mAodt^x%^ ^ -qA2_^qB X 



BO. 



(5.42) 



The right side of ( 5.42 ) contains the first time irreversible term in our conditions of motion. Its form and sign is a 
direct consequence of having solved our wave equations as an initial value problems, leading to the retarded solution 
in the asymptotic future. If they had been solved as a final value problem so that the solution would have been an 
advanced wave in the asymptotic past then one would have obtained a result similar to ( ^.42 ) but with the sign o f th e 
radiation reaction term reversed, thus maintaining the complete symmetry of the basic field equations (3.1) and ( [3.4| ). 
Put another way, the arrow of time points in both directions. Of course, if one were to replace t by —t in the past 
the advanced solution would transform into a retarded solution and the sign of the reaction term would be reversed. 



VI. ALD AND EIH COMPARED 



We have now arrived, starting from the empty-space Einstein-Maxwell equations (3.1), (3.4) and ( |3.5| ) and making 
only the ass umpt io n of s pheri cal sy mmetry for the sources of the gravitational and electromagnetic fields, at a set of 
constraints (5.15), (5.27) and ( 5. 42] ) on the motion of these sources. They are applicable only to slowly moving sources 
and then only in the asymptotic future of some initial value surface. As such, they stand in stark contrast to the ALD 
equations of motion (2.18) which arc held to be exact and to hold for all times. If one were to forget the problems 
associated with their derivation one would surely favor them over the EIH equations even on esthetic grounds. Except 
for one problem. The ALD equations admit physically unacceptable solutions - the so-called run-away motions. The 
EIH conditions on the other hand do not suffer from this defect as we shall see. 

The source of the run-away solutions to the ALD equations is of course the third derivative terms in the expression 
for the radiation reaction force appearing in (2.18). Third derivatives also appear in the radiation reaction term in 
( ^.42 ). However there is an essential difference in the two cases. In (5.42) the third derivative is with respect to 
ti and not to the single time variable t appearing in ( 2.18| ). The dependence of the source coordinates on ti (and 
ts) is fixed by the Coulombic and post-Coulombic conditions ( ^.15| ) and (5.27). Equation ( 5.42| ) serves to determine 
the dependence of the source coordinates on damping time and contains only first derivatives with respect to this 
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variable. As a consequence, the motions allowed by the EIH conditions are all well behaved. That there are no 
run-away motions allowed in general relativity follows from the work of Gibbon, Hawking, Horowitz and Perry |2^ ] 
who have shown that the Bondi mass of a system of charged masses and/or charged black holes is positive definite 
and the fact that if the Bondi mass is initially finite, its time derivative is negative semi-definite. The Bondi mass 
thus appears to play the role of a Liapunov function and its existence therefore strongly suggests that these systems 
are asymptotically stable. 

To better understand the fundamental difference between the ALD and EIH equations let us apply them to a simple 
problem - a non-relativistic oscillator. In this case the ALD equation reduces, for an electron, to 

X +ujIx — T X (6.1) 

where r = le^/mc"^ is the light travel time across the classical radius of the electron. In keeping with the assumption 
of non-relativistic motion I assume that ojqt = e <C 1. Finally, introduce a new independent variable i' = u>ot to 
obtain the equation 

x +x — ex. (6.2) 
The solution will be proportional to an exponential e^* where A satisfies 

+ 1 - eA^ = 0. (6.3) 
Rather than write out the exact solution for A one can construct a power series solution in e. One finds that 

A±=±z-ie + 0(e2) (6.4) 



There is a third solution given by 



\ = -+Oie) (6.5) 
e 



which is the run-away mode. It has the peculiar feature that A — s- oo as r — > 0, clearly indicating that something is 
wrong with our equation. Dirac sought to remedy the situation by requiring that the initial acceleration, which must 
be specified along with the initial position and velocity since we are dealing with a third order equation, should be 



such as to rule out the run-away solutions. |26| Aside from its ad hoc character, this requirement was shown to lead 
to acausal behavior. In the present case it would amount to taking the amplitude of the run-away mode to be zero. 

Consider now the solution of the same problem using the EIH conditions of motion. In this case ( p3.15 ) would have 
the form 

m2 x +LLi'g'^x = (6.6) 

where dots now indicate differentiation with respect to ti and a prime is used to differentiate EIH quantities from 
their corresponding ALD counterparts. Motions allowed by these conditions are of the form 

X — A±{t4)exp{±iu}'Qti) (6.7) 



where the amplitudes A± are now taken to be functions of t^. When this result is substituted into ( 5.42| ) we obtain 
the equation 

2iu;'o^dt,A± = -t'u;'o^A± (6.8) 

where r' is the same function of the scaled charge and mass of the electron as r is of their unsealed counterparts. 
With the scalings introduced at the beginning of Section IV we have e'^r' = r. Taking c'cu'q = cuq we find that 

A± = Ao±exp(i™^2t). (6.9) 



Thus we obtain two damped modes c orre sponding to the two roots (6.4) of (6.3). But we do not obtain a motion 



corresponding to the run-away mode (6.5) 
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VII. SUMMARY 



The conditions of motion derived above together constitute asymptotic conditions on the motion of charged, spher- 
ically symmetric, compact, slowly moving sources accurate to 0{S'^,e'^). Given the current positions and velocities of 
a collection of such sources sufficiently far into the future of an initial startup time they will predict their motion for at 
least a finite time into the future to this order of accuracy. They were derived from the source-free Einstein-Maxwell 
field equations without additional assumptions concerning force laws or the form of inertial forces. Their derivation 
did not require assumptions concerning which parts of the fields act on which sources as does the derivation given 
by Landau and Lifshitz, for example, nor did it require the renormalization of divergent integrals such as in the 
Dirac derivation. And these conditions do not admit non-physical solutions which violate the approximations used 
in their derivation. In this connection. Landau and Lifshitz are, to the best of my knowledge, the only authors who 
suggest that the existence of run-away solutions is evidence for the limited applicability of the radiation reaction 
force. But then, like other authors, they inexplicably ascribe this state of affairs to the renormalization of the infinite 
electromagnetic self energy of point sources when in fact no such renormalization is necessary. 

We see thus that the outstanding problems of classical electrodynamics, the need for infinite renormalizations and 
the run-away solutions, find their resolution within the framework of the general theory of relativity. In addition we 
see how irreversible effects, here radiation damping, arise from reversible fundamental equations. Rather than having 
to assume an outgoing wave solution in the slow zone valid for all times it is sufficient to assume only that the initial 
conditions are physically reasonable, that is, that the initial wave energy is finite. 

One can ask: is it possible to do better than what has been accomplished here? One can of course go to higher 
order in e in the 0{S^) approximation, obtaining thereby higher order special relativistic corrections. One can also 
go to higher order in 6 to obtain corrections describing electro-gravitational effects. But can one derive equations of 
motion for these sources rather than conditions of motion? And can one at least derive conditions of motion for 
arbitrarily moving sources that are special relativistic invariant. The answer to all of these questions is, I believe, 
most likely no. In the fast zone the motion of the sources will depend in a complicated way not only on their initial 
positions and velocities but on the initial values of the fields as well. Since during this period there will not in general 
be any small parameters associated with the motion one would have to solve the field equations exactly in order to 
evaluate the EIH surface integrals and that, so far, has not been done. It is only if one waits long enough for all of 
the initial transients in the motion to die down that one can expect to derive relatively simple conditions of motions. 
For a similar reason one cannot derive special relativistic equations or conditions of motion. All attempts to do so 
in the past have been singularly unsuccessful, both literally and figuratively. One of course should not feel to badly 
about this state of affairs since there really does not exist a self-consistent special relativistic description of interacting 
charged bodies. 

To conclude, I want to argue that the derivations of EIH and those here shed an entirely new light on the meaning 
of equations of motion which lie at the very heart of classical physics. It is not just the fact that conditions of motion 
can be derived from the empty-space field equations of general relativity that is important but even more so is their 
very nature. In contrast to equations of motion they are both approximate and asymptotic, they do not hold for all 
times and they cannot be used to solve initial value problems. This of course does not mean that everything done in 
the past in classical mechanics is wrong because it used these conditions of motion as equations of motion since their 
application has always been compatible with the approximations made here. But it does mean, I would also argue, 
that, in a very real sense, all of classical mechanics and electrodynamics can be looked on as a verification of general 
relativity since one does not have to postulate the conditions of motion which underlie these theories but can derive 
them from general relativity. 
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APPENDIX: 



In this appendix I w ill ou tline the method of determining the field components needed to derive the post-Coulombian 
conditions of motion (5.27). The first component needed is /i4°, which satisfies 



2, 00 



--ho 



(Al) 
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This equation can be solved exactly by 



1 -^^^t; ■ r • s -^^j. ■■ r 1 . '^A2 

where 



lY^TnA-A-y/ +^1 i;]+45]^ (A2) 

^ TA ^ TA 



YF = n\n% - and = n^, , n^, ^ (A3) 

3 TA 

are spherical harmonics of orders two and one, respectively. The coefficient mA4 in the homogeneous term in the 



solution (A2) is fixed by ( |3.q ) with = using the ©(e"") terms. The surface independent terms can be made to 
vanish in this order if we take dt^niAo = 0, set 



1 .2 , 1 ' 

= 7;^A2 Xao 

' ' ' XABQ 



and include the constant of integration in a redefinition of mA2- 

The next quantity we need to determine is ft,™". Unfortunately, there does not exist a closed form solution of 
equation ( 5.26| ) for these field components. However, all we need to know are those parts of /i™" required for the 



evaluation of the surface independent terms of the first integral in equation (pM) with n — s over a sphere surrounding 



the ^th source. These parts can be obtained by expanding the right hand side of equation (5.26) in a Laurent series 
in VA and keeping those terms that generate these desired parts. One finds in this way that 

h\'A~2mA2{5'''n\ Xaq-u^ i^o ""a + — (A5) 

TA 



where the ~ indicates equality modulo terms that do not contribute surface independent terms. The coefficient h 



ran 



in the homogeneous term is fixed by the imposition of the deDonder conditions (5.7) and is given by 

= 4m^2 a^Aoa^lo (A6) 
In a similar manner, the needed parts of h™^ are obtained from equation ( ^.25 ) and are found to be given by 



ft°l ^ 2qAn\ F12 + 2<Za(^;o "A" O F^^- 
-2mAQ — {n\n\ - ^5") x\^x\^x\ -2qAn\ (2 FJ^^-V i^2A)+ 

TA 3 

+ ^ (A7) 

rA 

The last term is a homogeneous term which must again be determined by the imposition of the deDonder coordinate 
conditions and which in turn requires the knowledge of the h^^ found above. When the deDonder conditions are 
applied one finds that 

Ja yn^Ao XaqXao ■ (A8) 
When these expressions for and are used t o eva luate U^'^ and the surface independent term s of the 



first integral in equation ( |3.6| ) together with equation ( |5.2l| ) one obtains the post-Coulombian conditions ( |5.27|) as 
described in the text. 
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